Geometric frustration inherent to the trilhum lattice, a sublattice of the B20 structure 
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We study the classical Heisenberg model on a recently identified three dimensional corner-shared 
equilateral triangular lattice, a magnetic sublattice to a large class of systems with the symmetry 
group P2i3. Since the degree of geometric frustration of the nearest neighbor antiferromagnetic 
model on this lattice lies on the border between the pyrochlore (not ordered) and hexagonal (ordered) 
lattices, it is non-trivial to predict its ground state. Using a classical rotor model, we find an ordered 
ground state with wavevector (^^,0,0) featuring 120° rotated spins on each triangle. However, a 
mean field approximation on this lattice fails to find an ordered ground state, finding instead a 
non-trivially degenerate ground state. As the mean field approach is known to agree with Monte 
Carlo on the pyrochlore lattice, the reasons for this discrepancy are discussed. We also discuss the 
possible relevance of our results to MnSi. 

PACS numbers: 75.10.Hk, 71.27.-fa, 75.25.-l-z 
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I. INTRODUCTION 

Systems with magnetic moments usually order as one 
approaches zero temperature. However, magnetic sys- 
tems that form triangle- or tetrahedron-based lattices are 
often able to avoid ordering until unusually low tempera- 
tures, and have been found to have interesting properties 
both theoretically and experimentallyiii^ii^iii^. Such lat- 
tices are generally referred to as geometrically frustrated 
as, for example, antiferromagnetic interactions between 
spins render the systems unable to find a unique ground 
state from magnetic considerations alone. It has been 
argued^ that the degree of frustration of a lattice can be 
quantified in terms of the connectivity of its triangular 
units. This leads to a hierarchy of lattices increasing in 
frustration from the ordered edge-shared hexagonal lat- 
tice, through the corner-shared triangle kagome, garnet 
and (3-Mn lattices, to the corner-shared tetrahedral lat- 
tice (common to the pyrochlore, cubic Laves and normal 
spinel structures). As we will show, the three dimensional 
corner-shared equilateral triangle structure of Fig. [1] is 
a fascinating new addition to this hierarchy. To avoid 
the cumbersome use of the terminology "three dimen- 
sional corner-shared triangle" (which also applies to the 
'distorted windmill' /3-Mn structure), we will henceforth 
name this lattice the trillium lattice^. The trillium lattice 
is a sublattice of the CO (B21), NH3 (Dl), NiSSb (FOi, 
UUmanite) and FeSi (B20) structures. For example, the 
Mn atoms of MnSi (B20 structure) form the trillium lat- 
tice. 

Recently Pfleiderer et al^ reported the observation of 
a "partially ordered" magnetic state in single crystals of 
MnSi near and above the critical pressure of this itinerant 
helimagnet. In this neutron scattering measurement, the 
well-defined magnetic Bragg peak {q ~0.0214 A~^) with 
wavevector along (111) was seen to shift to point predom- 
inantly along (110) (with q ^0.0304A~-'^) upon the appli- 
cation of pressure. The magnetic scattering pattern was 
seen to feature a sphere at fixed wavevector magnitude, 
with anisotropic intensities which weakly peaked along 
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FIG. 1: (Color online) The trillium lattice. Each site belongs 
to 3 equilateral corner-sharing triangles, {q, /3, 7, S} denote 
the four atoms of the cubic unit cell. Normal vectors to each 
triangle point along local (111) axes. (1,1,1) is a C3 rotation 
axis of the crystal; the front right lower triangle is related to 
the front left top triangle by translation. In blue (at left) we 
identify the layer number that each point of the equilateral 
triangles shown belongs to. The separation between layers 
4 and 5 is iuao/VS; between 3 and 4: (1-4m)oo/'\/3 and al- 
ternates. After triangles, the next smallest nearest neighbor 
connection has 5 sides as shown by the dashed line. 



(110). This sphere was found to onset at temperatures 
consistent with a continuation of the slope of T/v (where 
helimagnetism sets in). Coincident with the appearance 
of this unusual magnetic signature an anomalous non- 
Fcrmi liquid resistivity p ex T^'^ has been observed and 
found to persist up to twice the critical pressure^. The 
lack of a rapid recovery of p cx in such an apparently 
clean system (po O.SSfxflcm) poses an important chal- 
lenge to our understanding of Hertz-Millis-Moriya spin- 
fluctuation theory near a ferromagnetic or helimagnetioiS 
quantum critical point. 

Magnetic susceptibility measurements have been found 
to exhibit high temperature Curie-Weiss tailsii with 
Curie- Weiss constants tuning continuously from weakly 
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positive temperatures^, consistent with the onset of heh- 
magnetic order at ambient pressures, to weakly negative 
offsets above the critical pressure. This may be indica- 
tive of a crossover from ferro- to antiferromagnetically 
dominated fluctuations, although debate exists in the lit- 
erature on this point^®. Therefore, we ask whether any 
link might possibly be made between an enhancement of 
the antiferromagnetic couplings and the observed quasi- 
degenerate nature of the spin-spin correlations in this 
material at high pressures. 

To begin to address this question, we need to under- 
stand the role of antiferromagnetic interactions on this 
lattice. Here we advance a microscopic description of the 
physics on the magnetic lattice of MnSi. It is impor- 
tant to notice that the Mn bonds of this lattice form a 
trillium latticeiS. (see Fig. [1]) whose degree of geometric 
frustration lies on the border between previously stud- 
ied frustrated systems and those known to order despite 
frustration inherent to their lattices. We first classify the 
degree of magnetic frustration of the classical Heisenberg 
antiferromagnet on this lattice, and then use mean field 
theory and a mapping to a classical rigid rotor model, to 
explore the corresponding low energy physics. 

In our model, we have chosen only to include the local- 
ized spin degrees of freedom. For the majority of mag- 
netic systems forming on this lattice this is clearly an 
oversimplification, and a more complete model for the 
monosilicides might be the Kondo lattice model. Within 
this model magnetic correlations which arise between 
neighboring localized spins are mediated by conduction 
electrons. The basic magnetic physics of such an RKKY 
interaction-'^^ is thought to be captured by the extended 
Heisenberg model. This paper represents the first step to- 
ward such a treatment which will be presented in the near 
future^^ and is expected to yield a rich phase diagram and 
potentially account for the small positive Curie- Weiss off- 
set seen in MnSi near its critical pressure. 

Mapping the antiferromagnetic Heisenberg model to a 
classical rotor model, we find an ordered ground state 
with wavevector (|^,0,0) featuring 120*^ rotated spins on 
each triangle. Monte Carlo calculations have been carried 
outi^ and support this conclusion. However, the com- 
monly used mean field approximation does not find an 
ordered ground state. Rather, mean field theory finds a 
macroscopically degenerate ground state with degenerate 
wavevectors forming a sphere-like surface in momentum- 
space, reminiscent of the partially ordered state of MnSi. 
Although the detailed anisotropy of the structure fac- 
tor is not reproduced, to our knowledge this is the first 
model to find such a partially ordered ground state. The 
disagreement between Monte Carlo and mean field ap- 
proximations is surprising in light of the excellent quali- 
tative agreement between these methods on the kagome 
and pyrochlore lattices, and the reason for this disagree- 
ment is discussed. The strength of magnetic order within 
the rotor model and its relationship to the mean field 
results is presented, laying the groundwork for future fi- 
nite temperature studies where Monte Carlo and mean 



field results find qualitative agreement to reasonably low 
temperatures^^. 

The outline of this paper follows. In Section II an 
overview of recent excitement in the study of geometri- 
cally frustrated materials allows us introduce the trillium 
lattice in its proper context. In Section HI we employ the 
standar d^°i^^ mean field approximation which implicitly 
relaxes the constraint on the spin = 4) at each 

site and calculate the neutron structure factor. In Sec- 
tion IV we numerically impose a hard spin constraint 
(sf — 1) at each site. We use a classical rotor model to 
characterize the low-lying excitations and ground state. 
Section V contains a discussion of the main results of this 
work. 



II. GEOMETRIC FRUSTRATION 



A. The trillium lattice 
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{S+y,5+y+i,a, a+y,'y, 7.4) 



TABLE I: The nearest neighbor sites in relation to the unit 
cell depicted in Fig. [1] The subscript denotes translation by 
one unit cell in the direction noted. 

The B20 crystal structure to which MnSi belongs con- 
tains two interpenetrating sublattices. For MnSi, one 
sublattice consists of the magnetic Mn atoms, while the 
other consists of non-magnetic Si atoms. Each sublat- 
tice forms an infinite three-dimensional lattice of corner- 
shared equilateral triangles, which we have named the 
trillium lattice. The trillium lattice forms a simple cu- 
bic lattice with a four site basis. Elements of this basis 
are listed in Table |T] and shown in Fig. [T] First, sec- 
ond, and third nearest neighbor bonds are all found-^ to 
form networks of independent, corner-shared equilateral 
triangles, such that even in systems with dominantly fer- 
romagnetic correlations between nearest neighbors, frus- 
tration may play a role. The symmetry properties of this 
lattice P2i3(r^) are known to include C3 axes through 
the center of each triangle and a screw axis, as discussed 
for example in Bradley and Cracknell^^. 

The history of the corner-shared tetrahedral lattice 
suggests that it is worth investigating systems other than 
the one of current interest to lay the basis for future work. 
It surely was not appreciated by Pauling that his famous 
estimate^'^ of the entropy of water ice (a strongly polar 
molecule) might lead to the explosion of interest in anti- 
ferromagnetically correlated Heisenberg spin systems on 
this lattice. Harris et found that large spin systems 
with ferromagnetic correlations could exhibit spin ice be- 
havior with its large low temperature entropy. Such a re- 
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alization so many years after Andersoo^^ and Villain's^^ 
pioneering works on antiferromagnetic spins on this lat- 
tice, must have come as a surprise to the early work- 
ers on magnetically frustrated lattices. As Harris et 
al explained, spin ice physics can arise due to crystal 
anisotropics which can pin large spins to lie along local 
Ising directions-symmetries of the crystal structure. 

While the focus of this paper is on the antiferromag- 
netic Heisenberg model on the trillium lattice, we briefly 
summarize systems that form the trillium lattice. The tril- 
lium lattice is a sublattice of the CO (B21), NH3 (Dl), 
NiSSb (FOi, UUmanite) and FeSi (B20) structures. Al- 
though both CO and NH3 (Dl) are non-magnetic, they 
have large dipolar electric fields. Magnetically interest- 
ing B20 structure systems include binaries of most of the 
transition metal elements with Si or Ge, although struc- 
tures have also been reported^! with Sn, Al, Ga, Hg, Mg 
and Be as the second member in addition to ternary com- 
pounds. Among these, the "Kondo insulator"^ FeSi has 
achieved prominence as an important constituent of the 
earth's mantle. Fci-^^Coj^Si behaves as a doped magnetic 
semiconductor and exhibits a large anomalous Hall effect 
which has been argued to have spintronics applications^ 
and MnSi is an enhanced mass metal close to critical- 
ity. The FOi structure is believed to include EuPtSi and 
EuPdSi, spin ^ systems exhibiting (ferromagnetic) Curie 
laws from 300K down to 5 K. Mossbauer spectroscopy has 
shown the europium atoms to be in a divalent (Eu^+) if 
state'^'^. If crystal fields are such that this large moment 
aligns along the local (111) axes, this could be a good can- 
didate for the realization of spin-ice physics on the tril- 
lium lattice as we will outline in future work^. In light 
of the large number of systems which have shown some 
evidence for strongly correlated physics and the prospect 
for future systems to study, it is surprising that to our 
knowledge, this is the first treatment in terms of local 
moment based magnetism on the trillium lattice. 



B. Itinerant frustrated magnetism 

Upon encountering unusual magnetic signatures on a 
geometrically frustrated lattice, even for metallic sys- 
tems, it is common to isolate the magnetic degrees of 
freedom by first treating a Heisenberg model. A number 
of frustrated systems show metallic behavior with rea- 
sonably enhanced effective masses. Among these, the 12 
magnetic sites of /3-Mn (A15 structure) have been shown 
to form the distorted windmill (corner-shared equilat- 
eral triangle) lattice, which has a P4i32 symmetry. The 
Heisenberg model on the distorted windmill lattice has 
been studied at a mean field level and shown to possess 
a macroscopic ground state degeneracy along the (111) 
axisS. The magnetic sites of both LiV204 (the only 
d-electron heavy fermion) and (Yo.97Sco.o3)Mn2 reside 
on corner-shared tetrahedral lattices. As detailed below, 
this lattice has been the subject of intense theoretical 
study. To date all such strongly correlated frustrated 



metals have shown evidence of strong nearest neighbor 
antiferromagnetic correlations. 

Theoretically, the classical antiferromagnetic Heisen- 
berg model for 0{Af) spins was first solved on the kagome 
lattice where the results were reported to be exact in the 
A/" — > 00 limit'''^ and recently this treatment was extended 
to the corner-shared tetrahedral lattice^''. Isakov et 
showed that the spin correlations could be understood 
in terms of a ground state constraint requiring the total 
spin on each tetrahedron to vanish. These results agree 
remarkably well with classical Monte Carlo simulations, 
in particular reproducing a "bow-tie" -like structure first 
found in results of Zinkin et alM-. Qualitatively similar 
structures are seen in quantum treatments of the Heisen- 
berg model on the corner-shared tetrahedral lattice;^ 
Experimental evidence for such quasi-degenerate bow- 
tie structures has been found in the neutron scattering 
measurements^ of itinerant (Yo.97Sco.o3)Mn2. More- 
over, the authors^S have pointed out that these results 
are additionally applicable to both spin ice (which has 
net ferromagnetic correlations) on the pyrochlore lattice 
and cubic (water) ice. In this paper, we study a Heisen- 
berg model on the trillium lattice to understand the geo- 
metric frustration inherent to this lattice and its possible 
relevance to MnSi. 



C. Quantifying frustration 

For a true energetic minimum, we know® that the 
Heisenberg model can be expressed as, 

H = J Si ■ Sj 

A i 

= —bj per spin. (1) 

Here taken as the sum over all nearest neigh- 

bor sites, such that on the lattice each bond is dou- 
ble counted, Sa is the total spin on each corner-shared 
triangular (or tetrahedral) unit of the lattice, and b 
is the number of connected units at each site as illus- 
trated in Fig. O Here we have used the constraint, 
I Sip = 1. However, note that the same minimal en- 
ergy state is reached if one relaxes this constraint to, 

|sf P + |sf P + Wi? + = 1 where a, /?, 7 and S 

might represent different sites in the unit cell. Minima 
(for J > 0) occur when \Sa\'^ — and have energy -bJ 
(per spin). Note the difference between the minimal en- 
ergy state on the trillium lattice (-3 J) and that of the 
corner-shared tetrahedral lattice (-2 J). 

It is instructive to classify the degree of frustration 
of the trillium lattice in relation to previously studied 
lattices (see Table HI]). Counting of the number of de- 
grees of freedom minus the constraints to maintain the 
ground-state condition®, |S'aP = 0, we find, the number 
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(a) (b) 




FIG. 2: The number of units, b, common to a given site on 
common lattices: (a) kagome: 6=2; (b) hexagonal: b—6; (c) 
corner-shared tetrahedra: 6=2; and (d) corner-shared triangle 
(trillium, distorted windmill): 6=3. 



of degrees of freedom of the ground state, Z?, to be given 
by, 

for I Sip — 1, and 

for |sf p + |sf p -I- |s7P + |sf p = 4. Here q is the number 
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D 
N 


q 
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ground state 


corner-shared 
tetrahedra 


1 


4 


3 


2 


degenerate 


kagome 





3 


3 


2 


degenerate^i^ 


garnet 





3 


3 


2 


degenerate^ 


distorted windmill 


-1 


3 


3 


3 


degenerate at m^ 


trillium 


-1 


3 


3 


3 


this work 


hexagonal 


-2 


3 


3 


6 


non-degenerate 



TABLE II: Degeneracy of the ground state of the classi- 
cal Heisenberg model for several common triangle/tetrahedra 
based lattices. From Eq. [2] and Eq. |3l a naive count of the 
number of degrees of freedom available to the ground state 
helps to classify the degree of frustration of a lattice. For 
the trillium lattice, the hard spin constraint (rotor model) 
yields D/N = —1, which is found to order; the soft spin con- 
straint (mean field) leads to D/N = —1/4 which seems to 
exhibit a partial ordering. Note that the Heisenberg model 
is believed to order by disorder on the kagome lattice-^'*^ so 
strictly speaking having a macroscopic number of zero energy 
spin structures is not sufficient to avoid order. The spin- 
ice model on the trillium lattice will necessarily break the 
|SaP = condition, so Eq. [2] and Eq. |3]do not apply to this 
model and a larger ground state degeneracy is possible. 



of spins per triangular/tetrahedral unit, N is the number 
of units, and A/" is the number of spin components {Af=l 
Ising and J\f=3 Heisenberg) . We see that upon strongly 
enforcing the constraint in Eq. [5] (rotor, Monte Carlo) the 
Heisenberg model is expected to order. The mean field 
approximation enforces the relaxed constraint in Eq. [3] 
to produce a less negative value of D. From now on we 
will refer to the constraint of Eq. [2] as the hard spin 
constraint, while that of Eq. [3] will be referred to as a 
soft spin constraint. 

Prom this nai've counting argument, we would ex- 
pect to find a degenerate ground state only when D > 
0. However, it is known that this counting argument 
can go wrong if the constraints on the spins are not 
independent!^, such that a macroscopic degeneracy of the 
ground state can be present as seen in Table |lll Note 
that despite this ground state degeneracy, the classical 
Heisenberg model on the kagome lattice is believed to 
order by disorder in the T — > limit 4^"^ 



III. A MEAN FIELD APPROXIMATION 



A. A remarkable degeneracy 



The nearest neighbor Heisenberg model can be ex- 
pressed as, 



H = jJ2s,-s,^-Y^ S.MfS.,, (4) 

{ij) 9 



Here J >0 corresponds to antiferromagnetic interactions, 
N is the number of unit cells in the lattice and Sq — 

^ s^, s^, s'^, Sq y The four sites of the unit cell are 

labelled by {(5, a, 7,/?} as shown in Fig. [T] and Table [J 
The second equation has been obtained by the Fourier 
transformation, = i s^e"'''» , where k is chosen 

from {S, a, 7, /?} and rj'"^ demarks both the unit cell and 
the position within the unit cell, with. 
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FIG. 3: Dispersion curves for J >0, ao = 1. Notice the 
degenerate minimum at energy -3 J, and in the lowest graph a 
band crossing. The total energy spectrum is invariant under 
independent operations of —> —qx, qy —qy and qz 



where u = 0.138 for MnSi. 

Within the mean field approximation, the energy dis- 
persion curves, X°'{~q) can be obtained by diagonaliza- 
tion of the matrix Mq. Here a = {1, ..,4} labels the four 
bands arising from the 4-site unit cell. In Fig. [31 energy 
dispersion curves are plotted for various directions of ~q . 
The ground state is given by the minimum of A°("g'). A 
remarkable feature of this spectrum is that the minimum 
energy solution is degenerate, u-independent, and occurs 
on a sphere-like surface (see Fig. 0]) demarked by. 



cos ( ) + cos — ) -I- cos ( , 



9 

4' 



(6) 



where we have restored ao, the length of the unit cell 
for clarity. Hence, within this approximation, at T=0, 
the classical solution of the nearest neighbor Heisenberg 
model on the trillium lattice has a degenerate (or par- 
tially ordered) ground state. This apparent degeneracy 
is a hallmark of frustrated magnetism. 

Note that, for each direction of a minimum of the 
dispersion occurs at a single wavevector magnitude, but 
there are infinite directions of 'q satisfying Eq. [6] with 
the same minimal energy state, — 3J. The partial or- 
der of the classical ground state observed here (with a fi- 
nite ~q minimum) is qualitatively different from the com- 



pletely degenerate ground state observed in the corner- 
shared tetrahedral, distorted windmill and kagome lat- 
tices, whose dispersions show at least one flat band for 
their ground state2ii22. This implies that the trillium lat- 
tice is less frustrated than both the corner-shared tetra- 
hedral lattice and the kagome lattice. This was expected 
from the naive counting of Table [Til as three triangles 
must meet at a point vs. two for the kagome lattice. 
However, the different behavior of the distorted windmill 
lattice was quite surprising to us. The trillium and dis- 
torted windmill lattices have remarkably different disper- 
sions given the qualitative similarity of their local struc- 
tures. 

Additionally, it is curious to note that along certain 
directions^ in momentum space there is a dispersionless 
band at finite energies. Peculiarly, a band crossing occurs 
where this can become the lowest energy state for a given 
wavevector. If a similar band structure arises within a 
strongly correlated nearest neighbor hopping model, ex- 
tensions of this work might favor the creation of heavy 
mass particles at chemical potentials close to this energy. 
Such a hybridization model is common to most heavy 
fermion descriptions and the fiat bands of the corner- 
shared tetrahedral lattice have been argued to lead to an 
explanation for the unusual physics of the only d-electron 
heavy fermion, LiV204^. 



B. Spin-spin correlation functions 

To calculate the spin-spin correlation functions and 
neutron scattering structure factor within this approx- 
imation, it is helpful to write the partition function. Fol- 
lowing Isakov et al'^'^, we introduce Af spin components 
and constrain each site to have unit spins. 



(8) 



= J VSVXe 



Note that in moving from Eq. [71 to Eq. [9l we have 
implicitly made the mean-field assumption Ai = A for 
the Lagrange multiplier. This is in fact an assumption 
inherent to the treatment of the corner-shared tetrahe- 
dral lattice by Isakov et al2S, which was found to agree 
remarkably well with Monte Carlo results'^''^ . As the four 
sites within the unit cell are in principle different*'*, this 



6 



is a four-fold relaxation of the constraints on the system, 
and implies that we are only enforcing the single-spin 
constraint on average within the unit cell. Additionally, 
we have introduced a chemical potential fi in such a way 
as to shift the value of the minimal energy to 0, and de- 

fincd 5^ = s^'", s"^'", s'^'", s^'"^ To find the saddle- 
point solutions, we need to minimize Eq. [9] with respect 
to A. Introducing = S^U and U'^S'^g = 5^, to 
diagonahze M^" = U~^M^^U - /xl, we obtain, 

^ ^^^n,.adet(f(J6^ + zA))AA 

= y"x>Ae^^^-^^e'^o."^'"(^('^''+'^^) (9) 
I 



where e° = A° -f 3 J is the energy of the a-th band. The 
constraint on A is found to be, 4iV — V„ „ ml , x - To 
calculate the spin-spin correlators, we introduce an auxil- 
iary field - ( , p^''^ ) coupling to the 

spins as, 
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(10) 



and calculate (sk^^sJ^!^) — ^-—-^^-^-p—\p=o- After some 

algebra and a Hubbard-Stratonovich transformation, one 
arrives at. 



Z ^ VX 



where ,^ = pl^Uia and p'^^ ^g = J2t'=i Ujpl,_g. 

The saddle-point solution is then 
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As [/ is generally not real, we note that C/^^ — 



the quasi-spheres given by Eq. [SI However, the geome- 
try of the lattice leads to an unequal distribution of this 
weight as a function of angle, and (see Fig. 2]) a dis- 
perse angle-dependence. The maximal intensity of the 
structure factor on this spheroid is found along (110) 
and vanishes due to symmetries of the lattice near (111) 
and (100). However, there is actually very little weight 
within the entire first reciprocal lattice zone of this non 
Bravais lattice, as for the corner shared tetrahedral lat- 
tice. The sphere-like shape (not shown) centered about 
(^,0,0) has a distinctly different weight distribution, 
} (12) tuning from minimal weight near (-^,0,0) to maximal 



(11) 



C. The static structure factor 

The static structure factor can now be found as^, 
S{li) (X ^(s^4e*'i-(--'-')) (X ^(.«^,.«.-,), (13) 

under the assumption that for classical spins all spin com- 
ponents contribute equally. 

As T ^0+, the small \~q\ structure factor of the near- 
est neighbor Heisenberg model has all of its weight on 



near (j|^,0,0) with the absolute magnitude of the weight 
1000-folS larger yet qualitatively similar temperature de- 
pendence. Translation by a reciprocal lattice vector of 
the cubic unit cell finds another spheroid, with vary- 
ing weight and distributions. One is reminded of the 
disperse neutron scattering patterns seen on the corner- 
shared tetrahedral lattice^Sii^. There, "bow-tie" struc- 
tures appear as a result of interference between the two 
length scales given by the unit cell and the structure 
within the unit cell. On the trillium lattice one again 
has two length scales which leads to interesting interfer- 
ence patterns, but at low temperatures, energetic con- 
cerns restrict the neutron scattering weight to lie only at 
wavevectors given by Eq. [6] The convolution of this en- 
ergetic structure with the underlying interference pattern 
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FIG. 4: (Color online) Within the mean field approximation, 
the ground state is found to be degenerate around the sphere- 
like surfaces corresponding to solutions of Eq. [6] Here we 
show (the Qz = half of) the degenerate spheroid centered at 
q = (0,0,0). When T ^0+ and J >0, as shown in Fig. \5\ 
one expects at low temperatures a surface of magnetic Bragg 
peaks to arise in the neutron scattering structure factor coin- 
cident with the degenerate spheroids. Within each Brillouin 
zone, the relative weight of each of the Bragg peaks is given 
by the geometrical factor of Eq. [ill Interference between the 
two length scales (cubic unit cell, nearest neighbor distance) 
leads to unusual patternings of intensity of the static structure 
factor, <S'(q), around each degenerate spheroid. The pattern 
covering this surface in the first Brillouin zone shows areas 
of high intensity (purple/dark) along (110), moderate (yel- 
low) around the surface and vanishes along (100) and (111) 
(red). However, the overall magnitude of the signal is approx- 
imately 1000 times smaller than that seen in MnSi, and this 
pattern does not hold about the (|^,|^,0) lattice Bragg point 
(although the signal would be substantial there). 
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FIG. 5: (Color online) Along the qqO direction, the temper- 
ature dependence of the longitudinal structure factor. (Note 
that unlike in MnSi, the area under the curves seems to be 
approximately conserved for us), (inset) The structure factor 
is calculated along the longitudinal direction parallel to qqO 
as shown. 



produces the varying weights of the structure factor on 
the sphere-like shapes around different reciprocal lattice 
vectors. 

As in the kagome and corner-shared tetrahedral lat- 
tices, our model (with relaxed constraint) has a macro- 
scopic ground state degeneracy. To calculate the struc- 
ture factor, as shown in Eq. [12] and 1131 we do a simple 
sum over all classical spin configuration states weighted 
according to their energies. As the energy gap along the 
surface of the sphere-like shape is zero, at sufficiently 
low but finite tempertures, the system will predominantly 
sample the different ground state spin configurations. In 
Fig. [4] we plot the maximal value of S{q) as T ^ 0+ for 
each wavevector direction within the first Brillouin zone 
to illustrate the idea of degenerate spheroids. 

For the trillium lattice, the second and higher bands 
are reasonably well-separated from the lowest energy 
state. At sufficiently low temperatures compared to J, 
the lowest band dominates the accessible spin structures. 
One expects a strong favoring of a particular wavevector 
magnitude for each direction in momentum space. This 
is seen as a function of temperature along (110) in Fig. 

El 

It is instructive to investigate the properties of the ex- 
treme limit as T^ 0+ when only the ground-state spin 
configurations should contribute appreciable weight to 
the correlation functions. In this limit, analytic calcu- 
lations can be done^*^, as the structure factor becomes 
proportional to. 
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i=l 



(14) 



We see that the first term returns 1 as long as the 
spin eigenfunctions are normalized. The rest are ex- 
pressible in terms of the relative spins on each site in 
momentum space. For example, a q=Q ferromagnetic 
state would have simply Un = U^^ = ^, leading to a 
factor of 4 in total. One can show in this limit that 
Eq. [141 vanishes everywhere along (111) by symmetry, 
although local (111) states about the reciprocal lattice 
point [j^j Ja^)^^) vanish. The detailed nature of 

the structure factor pattern on the degenerate spheroid 
centered at (-^j 0) does not correspond to that seen 
in MnSi. 



D. "Spin" structures: the soft spin constraint 

It is useful to have a physical picture of how the clas- 
sical spins are arranged along the lattice in the degener- 
ate minimal energy ground state. From studies on the 
corner-shared tetrahedral lattice, we are used to the con- 
straint that all spins on each tetrahedron must add to a 
spin-0 state to be included in the degenerate ground state 
manifold. The analogous condition for the trillium lat- 
tice (as presented in Sect. II C) is for the spins on each 
triangle to add to zero. If we have a truly degenerate 
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ground state then at any finite temperature the classical 
spins will explore all these states. 

The eigenvectors of Eq. [4] hold information about the 
relative angles between spins. At any point in momen- 
tum space, if the eigenvector is non-degenerate then only 
one relative angle between spins is contained in this for- 
malism. Therefore, if a spin structure satisfying the hard 
spin constraint can be formed in the ground state, its 
spins will necessarily be coplanar. However, if such a spin 
structure cannot be formed in the ground state, then the 
degeneracy will be seen as an artifact of the soft con- 
straint approximation. 

To obtain a particular spin configuration of the ground 
state, we need to_perform a restricted Fourier transform, 

= jf ' ■ • s^|>. The relative position vectors are 

given by, rf = rf + (5, ^ - 2u, 1 - 2u), r] = rf + {1 ~ 
2m, i, i - 2u), and rf = rf + (i - 2m, 1 - 2u, i). For 
simplicity, we consider "spin" structures corresponding 
to the first Brillouin zone in just three high symmetry 
directions, those with q along (111), (110) and (100). 



1. (Ill) 




FIG. 6: (Color online) Looking down the (111) axis, "spin" 
configurations with the relaxed constraint of Sect. Ill A. 
Numbers denote the layer number (or height) within the 
structure, (a) The lowest energy qo = state along (111) is 
shown here. Dashed and dotted lines label triangles out of the 
plane. On these triangles the ground state has two antipar- 
allel spins, (b) A "spin" configuration at higher energy along 
(111) at the same wavevector, which coincidentally would be 
in the ground state manifold for either second or third near- 
est neighbor models. Note that one no longer has antiparallel 
spins on the small triangles of (a) out of the plane, but that on 
the dotted lines in (b) (corresponding to 3rd nearest neighbor 
bonds) one does have antiparallel spins out of the plane. 



To draw the corresponding "spin" structure with 
q along (111) at go = gf^i we take the hermitian 
conjugate of the right ground state eigenvector of Eq. 
[S] (which corresponds to S-q-q-q) to obtain, S-^ — 

(0' 73 ' ^^-^^^ - -here ^ = (3^^, gf^, 3^), 
which yields,{sf , sf , sj, gf }=-l=e^(f + 

{0, l,e^,e^}. The "spin" angle gains an additional 
factor of by the translation of one lattice constant 
along any of the {x, y, z} directions of the cubic unit cell. 
By inspection, we see that spin contributions on the a, /3 
and 7 sites are of equal magnitude, but are rotated by 
120" relative to one another. A pictorial representation 
of this state is shown in Fig. [H^a). 

In this structure every triangle is found to have total 
spin 0, as expected from Eq. [TJ However, the S site (la- 
belled 1,3,5 in Fig. [6] (a)) has spin 0. Triangles involving 
this site have their other two spins arranged in an an- 
tiparallel fashion. As discussed in Section IV D, unlike 
on the corner-shared tetrahedral lattice, on the trillium 
lattice one cannot add perpendicular spin components 
with a different q to restore the hard spin constraint 
while remaining in the ground state. 



2. (110) 

Within the mean field approximation, the 
largest neutron scattering structure factor in 
the first Brillouin zone was found to lie along 
(110). The minimal energy spin structure is found 

at qo where we find S-^ = 



-i^(l,2e*2ug+«7r^2e'(S"+5)'3-^^^e*(4«-i)9)^(-^ - ^) 

with -qo = (^cos-i(y^), ^ 008-1(^1), 0). We can 
express the "spins" in real space as {sf,s°',s],s'^} 

gi2cos ^(vT)}. A phase change of e^'^'^°^ ^(vl") (or a 
spin rotation of 75.5°) is gained by each site following a 
translation by one lattice constant in the x or y direction 
of the cubic unit cell. No change in the spin structure 
accompanies a translation by one lattice constant in the 
z-direction. This spin structure is plotted in Fig. [Tl^a). 

In this structure, the total spin on each triangle is 0, 
again satisfying Eq. [S] However, "spins" on the S and 
(3 sites are seen to have only half the magnitude of those 
on the a and 7 sites (refer to Fig. [1] for a definition of 
the {6, a, 7, (3} sites). 

3. (100) 

Along the (100) direction we find the energy min- 
imum to occur at on = — cos^^(i) = where 

aa ^2' 3ao ' 

- (i^-^^^^^^-^^^^^Wt- 9^). Here 
gjj — (^,0,0), implying that the real-space spin struc- 
ture is {4,sf,s7,sf} = ie'^'-" {l,e^,l,e*^}. A 
phase change of e*"f (or 120°) is gained by each site fol- 
lowing a translation by one lattice constant along the x 
direction of the cubic unit cell. Translations by one lat- 
tice constant along the y and z directions leave the spin 
structure unchanged. This spin structure is shown in Fig. 
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FIG. 7: (Color online) Looking down the (111) axis of the 
trillium lattice, two degenerate "spin" configurations of the 
relaxed constraint. Numbers denote the layer number (or 
height) within the structure. Dashed and dotted lines label 
triangles out of the plane, (a) The qqO state, (b) The gOO 
state. Note that the latter both forms a 120° state on each 
triangle and satisfies the hard spin constraint. 

WIM)- This is a spin structure with 120° rotated spins of 
unit magnitude on every triangle. 

IV. MAPPING TO A RIGID ROTOR MODEL 

A. The hard spin contraint 

On frustrated lattices, subtle effects arise from con- 
straints on the spins. On both the kagome and corner- 
shared tetrahedral lattices the gournd state is known to 
be degenerate within both the {JV oo) mean field 
approximation and {Af = 3) classical Monte Carlo. It 
is important to ask whether the ground state obtained 
by the mean field approximation on the trillium lat- 
tice can be reproduced by a method which satisfies 
the hard spin constraint. Within the mean field ap- 
proximation, the spin structures corresponding to the 
dispersion curves of Fig. [3] satisfy the soft spin con- 
straint. We have seen from Section III D, that "spin" 
structures pertaining to a single ordering wavevector, 
~q , within the mean field, feature coplanar spins, and 
that the magnitude of these spins is generically found 
to vary from site to site within the unit cell. However, 
the (^,0,0) spin structure of Fig. [7] (b) was found 
to satisfy the hard spin constraint. We see from the 
mean field that this structure, (and its equivalents: {{- 
|i^,0,0},{0,^,0},{0,-|^,0},{0,0,|2L},{0,0,-^}}), is 
the only coplanar member of the ground state of the clas- 
sical Heisenberg model satisfying a hard spin constaint. 
To fully investigate the low energy properties of the 
Heisenberg model while strictly enforcing the hard spin 
constraint at each site, it is useful to consider the spins 
to be rigid rotors. This allows us to numerically both: 
(a) investigate finite size lattices with periodic boundary 
conditions to determine whether the (100) ground state is 
truly unique or if non-coplanar spin structures (multiple 
wavevector) are able to satisfy the hard spin constraint 




FIG. 8: (a) To find the energy of a spin structure which is 
periodic along a particular direction, we design a finite size 
system satisfying periodic boundary conditions. For a finite 
size spin lattice with {L^, Ly, L^) unit cells in the {x,y,z) 
directions, where LiUo is the length of the system in the i 
direction, the spin structure repeats after {Lx/n, Ly/m, Lz/l) 
unit cells where {n, m, 1} are integers greater than or equal 
to one. Wavevectors with this periodic boundary condition 
satisfy q = (^2I2_ 27rm_ 22lL). We then minimize the energy 
of randomly oriented spins on this lattice which chooses one 
(or more) value for each of n, m and I. As an example, the cu- 
bic lattice drawn above repeats after (L, 2L, 1) unit cells and 
allows wavevectors of the form, q = (|^,|ff^,0). (b) To 
only consider spin structures which repeat the same number 
of times in each (non-trivial) lattice direction {n = m — l),we 
can write a one dimensional lattice with all unique unit cell 
positions and connections. For example, for our cubic lattice 
which repeats after (L, 2L, 1) unit cells, we can ask that the 
spin structure after a translation by one unit cell along the x 
direction is the same as the spin structure after a translation 
by two unit cells in the y direction. Wavevectors, q, with this 
additional periodic boundary condition satisfy n = m as de- 
fined above, and correspond to a particular wavevector direc- 
tion (in the example given here this is (g, |, 0)). Minimization 
of such spin structures allows us to calculate the dispersion of 
the model with a hard spin constraint for several directions 
of ordering wavevectors. 

while remaining in the ground state; and (b) investigate 
the low-lying energetic saddle-points for various wavevec- 
tor directions. 



B. Method 

One way to strictly enforce the constraint of unit spins 
at each site is to represent the spins in terms of rigid 
rotors with fixed magnitude and two angles, 9 and 0, 
where ~Si—{sin{6i)cos{<j)i), sin{9i)sm{(j)i),cos{6i)). We 
minimize Eq. [3] on finite size clusters in real space. With 
N spins in a finite size lattice, the link matrix correspond- 
ing to the nearest neighbor couplings of the Heisenberg 
model is an A'^ X matrix. Such matrices can be simply 
written from the connections of Table [J While this leads 
to large and rather unwieldy matrices, one can use pe- 
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riodic boundary conditions to reduce the large matrices 
to a smaller one-dimensional form, as shown in Fig. [S] 
Varying the number of unit cells, L, in a manner con- 
sistent with the periodic boundary conditions allows one 
access to a finite but non-ncgiigiblc number of low en- 
ergy points in momentum space. Accessible states have 
components of the form qi = -jj^- For each wavevector 
direction, q, we minimize the energy with respect to the 
magnitude of "g*. We have chosen 7 unique directions for 
comparison with the mean field results. 

C. A unique (^,0,0) ground state 

Within the rotor model we have found that the ground 
state of the antiferromagnetic Heisenberg model on the 
trillium lattice is ordered along (100)^ with wavevector 
go = ^ and features 120" rotated spins on every tri- 
ang le. "the energy of this state is found to be — 3J as 
expected from Eq. [1] The rotor model has the advan- 
tage that it gives us direct access to the spin structure, 
which is found to have the same form as shown in Fig. [7] 
(b). We note that this is the spin structure of the only 
member of the mean field ground state which satisfies the 
hard spin constraint. 

1. Energetic saddle-points: a sphere-like surface 

Any spin configuration within the rotor model auto- 
matically satisfies the hard spin constraint, so it is in- 
teresting to ask which features of the mean field approx- 
imation survive. From the mean field, we found that 
restricting the angles to lie within a solid angle bounded 
by (100), (110) and (111), allowed access to all distinct 
energetics associated with the ground state. It is thus 
natural to tile this area and ask: (i) whether for each 
wavevector direction the energetic saddle-point occurs at 
approximately the same wavevector as given in EqlHl (ii) 
about the relative magnitudes of the low lying excitation 
energies along the longitudinal (parallel to ~q) and trans- 
verse (perpendicular to ~q ) directions; and (iii) whether 
there exists a simple analytic relationship between the 
mean field and lowest energy real spin structures. 

Table IIIII presents the results extracted in this fash- 
ion for 7 directions, (i) We see that the ground state of 
—3 J is realized along (100) and that the extracted values 
of |~g*|TOior closely follow those of the mean field, form- 
ing a sphere-like surface of saddle-points. (We find that 
along (100) the spins are always coplanar so the mean 
field values are correct), (ii) Remembering that the point 
(^,0,0) has the energy -2-^/2 J, we see that spin excita- 
tions along (100) are quite soft and of a similar order of 
magnitude to those seen around the sphere-like surface 
of saddle-points, (iii) In the Section IV D we will present 
an analytic derivation of the minimal energy state along 
the qqq symmetry axis of the crystal which corresponds 
to adding a ~q=Q state to the mean field result. 



dir 


Eniin 


\nr' 

\ y 1 rotor 


— ^(2) 
Q rotor 


l'?|mf 


<im 


-2.6786J 


^ - 1-814 


(0,0,0) 


TV 


qqO 


-2.6666J 


1.7198 


? 


1.864 


qOO 


-3J 


^ « 2.09 


none 


2n 
3 


y 2 2 


-2.7632J 


1.776 


none 


1.862 




-2.874J 


1.894 


(47r, 27r, vr) 


1.902 




-2.9302J 


1.916 


(27r,7r,0) 


1.926 


111 


-2.8525J 


1.767 


(27r,27r,7r) 


1.835 



TABLE III: We see that the actual energetic minima along 
several wavevector directions continue to lie close to their 
mean field values, delineating a spheroid of saddle-points con- 
nected to the true mininum at ;^00. For simplicity we have 
set ao = 1. 



2. (100) saddle-points: a cross-check 



lattice dimensions 


T^rotor 






2x2x2 


-2.8284271 


-2\/2 


n 
an 


1x4x4 


-2.9449468 


-2.9449473 


TV 

2a„ 


4x5x5 


-2.9691749 


-2.969174998 


47r 
5an 


1 X 10 X 10 


-2.9915323 


-2.9915325 


Sir 
San 



TABLE IV: Periodic clusters of up to 400 rigid three di- 
mensional rotors of various shapes and sizes have been nu- 
merically minimized (see Fig. |8] (a)) for the antiferromag- 
netic Heisenberg model. For representative clusters, we show: 
Lx X Ly X Lz, the number of 4-site unit cells considered along 
each axis {x,y,z} of the cubic lattice; E^°l°'^ , the minimized 
rotor model energy; the lowest energy (100) state consistent 
with the boundary conditions; and Qmin, the magnitude of 
the wavevector of this (100) state. In all cases the mini- 
mized energy corresponds to the lowest energy state which 
orders along (100) consistent with the boundary conditions. 
Note that, in the last case, the boundary conditions allow 
the states with q = (0, which satisfy Eq. [6] and 

were found to belong to the ground state in the mean field 
approximation, but the spin structure minimizes to one with 
q = (0, 0, This shows that, with a hard spin constraint, 

low- lying excited states along (100) remain lower in energy 
than the saddle-points of the no longer degenerate spheroids. 

To check that the observation of order along (100) 
within the rotor model is not an artifact of the map- 
ping procedure shown in Fig. [8] (b), we have also con- 
sidered finite size lattices as shown in Fig. [S] (a). The 
minimization of such lattices allows a spin structure to 
explore various multiple wavevector states that might be 
missed by the restricted geometry of Fig. [8] (b). In Ta- 
ble [IV] we compare the energetics from minimizations of 
the three dimensional spin structures of various shapes 
and sizes. An analytic expression for the dispersion along 
(100) which satisfies the hard spin constraint is given in 
Table |Vl This has been used to calculate the third col- 
umn of Table IIVI at the wavevector magnitude denoted 
Qmin- For all studied cases, the lowest energy spin struc- 
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ture of the rotor model is consistent with the lowest en- 
ergy (100) state permitted by the boundary conditions. 
This shows that not only is the (^,0,0) ground state 
unique, but the longitudinal spin excitations along (100) 
are quite soft and dominate the low energy dispersion. 



a 


Eigenvalue (A*",^) 


Eigenvector (V'^oo) 


1 


-2J(cos(f)+2sin(|)) 


i(l,-l,ie-'(2u-i)g^_jg->(2„~l)9) 


2 


2J(cos(f)-2cos(f)) 


-(1, 1, — e~''^"~3'9, _e~*(2"~i)9) 


3 


-2J(cos(f)-2sin(f)) 


i(l,-l,-je-"(2«-i)<7^je->(2«-i)9) 


4 


2J(cos(2) + 2cos(2)) 


i(l,l,e-'(2"-3)^e-'(2"-4)9) 



TABLE V: Analytic eigenvalues and eigenvectors of Eq. |4] 
along (100). These eigenvalues satisfy the hard spin constraint 
at every site and the lowest energy spin structures correspond- 
ing to the eigenvectors directly correspond to those found in 
the rotor model. The longitudinal spin excitations from the 
ground state along this direction are at least as soft as exci- 
tations transverse to {j^, 0, 0) as explained in Table HVl 



D. Analytics: mean field to rotors 

Since real spins live in three-dimensional space, spec- 
ifying one angle in principle only tells us information 
about the relative orientation of two components of the 
spin. We see from Fig. [H] (a) that, with a soft spin con- 
straint, the mean field ground state with (111) order had 
120° rotated spins on triangles in the plane with nor- 
mal vector (111) and no spin at the 6 sites. All other 
triangles featured two opposite spins. It is instructive 
to add spin components normal to this plane to recover 
equal spins at each site. That is, to obtain a spin struc- 
ture with unit magnitude spins at each site, we perform 
a vector summation of the mean field spin structure of 
Fig. [6] (a) and spin components orthogonal to this plane. 
Clearly one must add a full spin to the S site, while con- 
tributions to the {a, 7, /3} sites must be added such that 
the spin components at each site square to one. On the 
{a, 7,/?} sites, we assign an amplitude sin(0) to spins of 
the mean field structure and an amplitude cos{9) to the 
spin component perpendicular to the plane. We refer to 
this perpendicular component of the spins as a g = 
component because, despite structure within the 4-site 
cubic unit cell, it is unchanged from one unit cell to an- 
other. This is shown pictorially in the inset to Fig. [5] on 
the four different magnetic triangles of the trillium lat- 
tice. Then we can write the total energy per unit spin as 
a function of the angle 9 as. 



-9sm^{e) +2(6cos2(6l) 4-6cos(6l)) 



(15) 



and to find a true minimum of this function, we simply 
solve for the derivative. 



de 



21sin(6') cos(6l) + 6 sin(6l) 




0.3 0.35 
q (in units of Ji/a) 



0.45 



FIG. 9: (inset) An illustration of the spin components of 
the noncoplanar spin structure at (gf-, gj-, 3J-) on the four 
unique triangles of the trillium lattice. Analytic results have 
been derived by the addition of a perpendicular q=0 spin 
component (see main text) to the mean field results of Fig. 
|6] (a), to enforce the hard spin constraint. In the main figure 
we present a comparison of the rotor model results (points) 
with analytic results of Eq. [T7] obtained by minimization of 
as explained in the text. The N x N link matrix written 
above each point shows the number of spin sites of each one 
dimensional lattice under consideration. The corresponding 
wavevector along (111) can be extracted as qi^Jtor ~ 
where n is an integer. 



(JV/4)ao ' 



which yields 9 = coa~^{=^). Substitution of 9 in Eq. [15] 
yields an energy of — 2.6786J. This is the same as the 
lowest energy solution found by the rotor model along 
(111), and is the true minimal energy state with this 
wavevector direction. 

Within the mean field approximation, the lowest en- 
ergy solution of Eq. d] can be shown to be Emf (qqq) = 

-2J(cos2(|) + ^3 cos2(§)sin^(f)). The spin structure 
remains as presented in Section III D, with the excep- 
tion of the relative angle change of the spin from one 
unit cell to the next; for generic values of q the triangles 
which include the 6 site no longer feature parallel spins. 
As we have analytic expressions for our eigenvalues, it is 
interesting to write the general expression in the (111) 
direction following the same steps, 



3(1 -( 



2(cos(§)cos(f + f )-l 



;f)cos(|)x 



1 



q 2tt 

cos(- + y) + 3( (^^) cos(f + f ) - 1) 



2(cos(§)cos(f + ^)-l)' 



(17) 







(16) 



which is plotted in Fig. [§] and agrees very well with the 
rotor model. 

Note that in principle this must be done also within the 
mean field solution of Reimers' paper—, in such a way as 
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to satisfy the hard spin constraint on the pyrochlore lat- 
tice. That is, an analogous mean field solution is found 
along (qqq) a.t q — in the pyrochlore case, but the high 
degeneracy of the lattice allows the addition of a q = 
(ie. the same in each unit cell) components to satisfy the 
hard spin constraint. Relative to mean field calculations, 
Monte Carlo calculations (which satisfy a hard spin con- 
straint) on the pyrochlore lattice would be expected to 
differ by a small factor. That this has not previously 
been appreciated is likely due to the high degeneracy of 
the pyrochlore lattice-it is not nearly as important an 
effect there as it is here. 



V. CONCLUDING REMARKS 

Motivated by recent neutron scattering and spin sus- 
ceptibility measurements of MnSi, we have considered the 
role of antiferromagnetic interactions and their possible 
relevance to a partially ordered state in this system. To 
address this issue, we have studied a Heisenberg model 
on the trillium lattice, a new three dimensional corner- 
shared triangle lattice formed by the Mn atoms of MnSi. 
The trillium lattice is also a sublattice of other systems 
including the CO (B21), NHg (Dl), NiSSb (FOi, UUman- 
ite) and FeSi (B20) structures. 

The Heisenberg model is an oversimplified model for an 
itinerant system such as MnSi. However, this model has 
been used to reproduce the measured spin correlations 
in other itinerant systems including (Yo.97Sco.o3)Mn2, 
LiV204 and /3-Mn, where it has been considered the first 
step in an understanding of the magnetic correlations 
present. We found that a certain degree of geometric 
frustration does exist on the trillium lattice with a naive 
estimation of the number of degrees of freedom available 
to the ground state less than the kagome and pyrochlore 
lattices, but on a par with the /?-Mn lattice and greater 
than the hexagonal lattice. 

We have treated the classical antiferromagnetic Heisen- 
berg model on the trillium lattice within the mean field 
approximation and by mapping the model to a lattice of 
rigid rotors. Within mean field theory on the trillium lat- 
tice, this model was found to have a degenerate ground 
state with wavevectors found to lie around a sphere-like 
surface given by Eq. [HI At small but finite temperatures, 
the neutron scattering structure factor would then be 
expected to feature weight predominantly confined to lie 
on this surface. Along this sphere-like surface, a disperse 
weight modulation would arise solely from geometrical 
factors and be roughly independent of temperature as 
found of the "bow-ties" of the corner-shared tetrahedral 
lattice. Measurements normal to these sphere-like shells 
would show resolution limited sharp peaks as the tem- 
perature decreased as shown in Fig. [31 These features 
bear a strong qualitative resemblance to neutron scatter- 
ing measurements near the critical pressure of the itin- 
erant helimagnet, MnSi, although notable differences are 
present. Chief among these are the radius of the observed 



sphere-like shape, which is smaller by a factor of 10, and 
the detailed nature of the pattern seen. More precisely, 
the measurements^ appear to have been made with re- 
spect to the lattice Bragg peaks at (f^,|^iO). About this 
point there is no reason to expect an inversion symmetry 
of the structure factor and the C3 symmetry of the lat- 
tice which is responsible for the vanishing of weight along 
(111) shown in Fig. |3]would not be present. 

By mapping the classical antiferromagnetic Heisenberg 
model to a rigid rotor description with a spin of unit mag- 
nitude at each site of the trillium lattice, we have found 
an ordered ground state at wavevectoi™ (-1^,0,0) fea- 
turing 120" rotated spins on each triangle. The rotor 
model is a real space construction which imposes peri- 
odic boundary conditions on finite size clusters of three 
dimensional vectors representing the spin at each site of 
the lattice. For several small cluster sizes (see Table IIVI) , 
we have found the minimal energy state to be consis- 
tent with order along (100) as analytically described in 
Table |Vl indicating that low-lying longitudinal fluctu- 
ations from the ground state are quite soft. We have 
designed an extension of the rotor model in such a way 
as to allow only one particular wavevector direction at 
a time. That is, rather than studying an M x x P 
lattice, where {M,N,P} are integers counting the num- 
ber of cubic unit cells in each direction, we have used 
symmetry to map the lattice to a one dimensional lat- 
tice of unit cells. For example, to investigate spiral order 
along (111), unit cells at {(1, 0, 0), (0, 1, 0), (0, 0, 1)} are 
all identical and can be replaced in the mapping by a 
single unit cell. Despite the strength of (100) order, this 
allows us access to low-lying energetic saddle-points. In 
Table IIIII we have shown that the wavevector associated 
to these saddle-points continues to lie close to solutions of 
Eq. [6l An analytic derivation of the wavevector and spin 
structure at the saddle-point along (111) has been devel- 
oped in Section HI D which illustrates the relationship 
between mean field and rotor model results. 

The discrepancy between the results obtained via the 
mean field approximation and those obtained by the 
mapping to the rotor model arises from the use of soft 
and hard spin constraints in the respective approxima- 
tions. The soft spin constraint allows the magnitude of 
the spin to vary within the unit cell as illustrated by the 
"spin" structures of Section HI D. For itinerant magnetic 
systems, it is not unusual for the magnetization density 
to vary spatially within an effective magnetic model. The 
hard spin constraint of the rotor model sets the magni- 
tude of the spin at each site to be 1. 

It is interesting to note that problems associated with 
the soft spin constraint of the mean field approximation 
arise in earlier treatments of frustrated lattices, yet do 
not seem to have been commented on. This is due, in 
part, to the extensive ground state degeneracy of these 
systems. For example, on the corner-shared tetrahe- 
dral lattice (a sublattice of the pyrochlore structure) , the 
mean field approximation yields^ a doubly degenerate 
fiat band with a three band degeneracy at q = 0. Even 
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though the coplanar mean field ground state spin struc- 
ture along (111) does not satisfy a unit spin magnitude 
(s^ = 1) at one of the four sites on each tetrahedron, 
it is possible to add spin components (with q = 0) to 
every site to recover this constraint. In this way, one 
finds a spin structure which belongs to the ground state 
and has components of its spins which rotate along (111). 
Calculations of the neutron scattering structure factor^° 
within the mean field approximation make the implicit 
assumption that the entire spin rotates. That along cer- 
tain wavevector directions only a component of the spin 
is rotating must imply that the detailed nature of the 
relative neutron scattering weight is not fully captured 
within the mean field approximation, despite the remark- 
able qualitative agreement seen with Monte Carlo results. 
Discrepancies between the soft and hard spin constraints 
are simply more clearly evident on the trillium lattice, as 
satisfying the hard spin constraint lifts the degeneracy of 
the ground state. 

The advantage of the mean field approximation is that 
it allows one to access all wavevectors and easily gen- 
eralizes to finite temperatures. Indeed, at temperatures 
T > 0.4 J one sees, from the rotor model (see Table HIH) . 
that all wavevector orientations should become accessible 
to the system. One might hope at such temperatures to 
recapture many of the qualitative features of the mean 
field description. A finite temperature comparison be- 
tween the mean field approximation and Monte Carlo 
results will be presented in the near future^^. Future 
worki^ will treat the extended Heisenberg model which 



is expected to exhibit a rich phase diagram, possibly rel- 
evant to currently available magnetic monosilicides. 

It is worthwhile to mention the possibility of frus- 
trated ferromagnetically coupled spin systems forming 
on this lattice. Indeed, on the pyrochlore lattice there 
are arguably more candidates for spin-ice physics than 
for degenerate antiferromagnetic ground states. As we 
will show in future work^, if Ising spins are restricted 
to lie along the local (111) axes of the crystal (this is 
the direction of the closest non-magnetic atoms in some 
of the UUmanites and may arise due to crystal/ligand 
field splittings), then antiferromagnetically coupled spins 
will be expected to order while ferromagnetically coupled 
spins will be extensively degenerate possessing a general- 
ized version of the two-spins out, two spins in rule of the 
corner-shared tetrahedral lattice. 
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